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ABSTRACT 

We present the discovery of a cold massive dust disk around the T Tauri 
star V1094 Sco in the Lupus molecular cloud from the 1.1 millimeter continuum 
observations with AzTEC on ASTE. A compact (r <320 AU) continuum emission 
coincides with the stellar position having a flux density of 272 mJy which is largest 
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among T Tauri stars in Lupus. We also present the detection of molecular gas 
associated with the star in the five-point observations in 12 CO J=3-2 and 13 CO 
J=3-2. Since our 12 CO and 13 CO observations did not show any signature of a 
large-scale outflow or a massive envelope, the compact dust emission is likely to 
come from a disk around the star. The observed SED of V1094 Sco shows no 
distinct turnover from near infrared to millimeter wavelengths, which can be well 
described by a flattened disk for the dust component, and no clear dip feature 
around 10 /im suggestive of absence of an inner hole in the disk. We fit a simple 
power-law disk model to the observed SED. The estimated disk mass ranges from 
0.03 to >0.12 M Q , which is one or two orders of magnitude larger than the median 
disk mass of T Tauri stars in Taurus. The resultant temperature is lower than 
that of a flared disk with the well-mixed dust in hydrostatic equilibrium, and is 
probably attributed to the flattened disk geometry for the dust which the central 
star can not illuminate efficiently. From these results together with the fact 
that there is no signature of an inner hole in the SED, we suggest that the dust 
grains in the disk around V1094 Sco sank into the midplane with grain growth 
by coalescence and is in the evolutional stage just prior to or at the formation of 
planetesimals. 

Subject headings: stars: circumstellar matter — stars: individual(V1094 Sco - 
solar system: formation 



1. Introduction 



Protoplanetary disk around pr e-main sequence sta rs are potential sites for planet for 



mation. In the standard model (e.g. lHayashi et al.lll985l ). after settling into equilibrium, the 



gas and solid components in the disk are separated and dust grains begin sinking to the mid- 
plane to form a thin dust layer which eventually fragments into planetesimals, ingredients of 
planets. The radial profile of the surface density structure at this stage is li kely one of the 
most important parameters to control what kinds of planets form in the disk (IKokubo &: Ida 
2002h . 



Since the gas and dust are well mixed in the beginning, the grain growth and sedimen- 
tation are key processes to separate dust from gas. Theoretical studies have suggested that 
dust grains are removed from the flaring surface and will sink to the m i dplane with growin g 
by coalescence along the disk evolution ( IDullemond fc Dominikl 120041 ; iTanaka et al.ll2005l ). 
One way to find such evolved disks is to carefully examine the Spectral Energy Distribution 
(SED) of young stars as well as detection of (sub) millimeter continuum emission, because 
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the shape of the SED of young stars is sensitive t o its temperature structur e mainly deter- 



mined by the vertical distribution of dust grains (jChiang fc Goldreichl 119971 ). iTanaka et al. 



( 120051 ) studied the dust growth and settling in passive disks to investigate the evolution of 
disk structure and SEDs. Their numerical results indicated that the flux densities at mid- 
infrared (MIR) from a disk get lower as a result of dust grain growth and settling from the 
surface layer. 

In the last two decades, extensive observations from optical to millimeter wavelengths 
have been conducted to study disk properties ultimately determining the nature of forming 
planets. Pre-main sequence stars, particularly in the quiescent phase such as weak line T 
Tauri stars (WTTSs) or Class III stars are good candidates for searching such evolved dusty 
disks in the process of sedimentation. A pioneer survey in (sub) milli meter continuum emis- 



sion of the circumstellar dust aroun d young stars was conducted by iBeckwith et al 



and followed by many other groups (lAndre fc Montmerle 



1994 



Qsterloh fc Beckwith 



Nuernberger et alil997l : lKitamura et al.ll2002HAndrews fc Williamsll2005l . l2007h . 



1990) 



1995; 



Kitamura et al 



(120021 ) made an imaging survey of single T Tauri stars in Taurus and derived the physical 
properties of the disks combined with the SEDs. They found the radial expansion of the 
disks attributed to outward transp ort of angular momentum w ith decreasing Ha line lumi- 
nosity, a measure of disk evolution. Andrews fc Williams! ( 120051 ) made a large survey of disk 
emission of more than 100 young stars in Taurus and found that there is a distinct difference 
in detection rates between classical T Tauri star (CTTS; 90%) and WTTS (15±5%). No 
massive evolved disk around WTTS, however, was found. 

In the present paper, we report the discovery of an unusual cold massive disk around 
the T Tauri star, V1094 Sco, at the southern tip of Lupus 3 at a distance of 150 pc. The 
source has a quiescent accretion disk presumably in t he transient stage from CTTS to WTTS 



becau se its equivalent width of Ha ranged 7-14 A ( Krautter et al.lll997t IWichmann et al. 



19991 ). V1094 Sco is known as a single star whose sp ectral type, mass, and age are K6, 0.78 
Mq, and 3 Myr, respectively ( IWichmann et al.l 119971 ). 



2. Observations 

Our observations have been performed with the Atacama Submillimeter Telescope Ex- 
periment (ASTE), the 10 m submillimeter telescope operated by Nobeyama Radio Observa- 
tory (NRO), a branch of National Astronomical Observatory Japan (NAOJ), in collaboration 
with University of Chile, and Japanese institutes including University of Tokyo, Nagoya Uni- 
versity, Osaka Prefecture University, Ibaraki University, and Hokkaido University. The tele- 
scope is located at Pampa la Bola, Chile, an altitude of 4860m, and is remotely operated from 
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the ASTE ope ration facilities in Chi le and Japan through the network observation system 
N-COSMOS3 (IKamazaki et al.ll2005l ). We em ployed two receiv ers, AzTEC 1.1 mm camera 
( Wilson et al. 2008 ) and CATS345 receiver ( Inoue et ali 2008 ). for millimeter continuum 
observations and molecular line observations at submillimeter wavelengths, respectively. 



2.1. Millimeter Continuum Observations with AzTEC 



AzTEC is a 144-eleme nt bolometer array camera installed on the ASTE between April 
2007 and December 2008 (IWilson et al.ll2008l ). In this period, the array was configured to 
operate at a wavelength of 1.1 mm. The angular resolution and field of view of AzTEC on 
ASTE were 28"and 7'8, respectively. 

We have made wide- field (40'x40') millimeter continuum observations of the Lupus 3 
molecular cloud including the target, V1094 Sco, in July 2007. The individual on-the-fly 
maps were obtained in the raster-scan mode with a scanning velocity of 100" s -1 . The 
telescope was scanned along the AZ-EL coordinates with a spacing of 1'95, a quarter of 
the field of view. All the data were taken at night or early morning and the average 270 
GHz radiometric opacity at zenith was ~ 0.05 during the observation period. Telescope 
pointing was checked by mapping the quasar J1517-243 every one and half, or two hours. 
The pointing corrections are applied during the data analysis. The flux conversion factor 
from optimal loading to source flux density of the AzTEC array was calibrated by observing 
Uranus, Neptune, or 3C279 every night. 



The AzTEC data set was reduced by using the AzT EC Data Re d uctio n Pipeline (IScott et al 



20081 ) for the identification of point sources. Following IScott et al.l ( 120081 ) . we employed the 
Principle Component Analysis (PCA) cleaning method to remove the atmospheric emission. 
The pixel size of the final map was set to be 3", 10 % of AzTEC beam size. The noise level 
and effective resolution of the final map are 5.0 mJy beam -1 and 35", respectively. 



2.2. Sub-Millimeter Molecular Line Observations with CATS345 

12 CO and 13 CO J=3-2 line observations were performed in March 2008 with a 2 sideband- 
separating (2SB) receiver, CATS345, on ASTE. The half power beam widths (HPBWs) at 
345.79 and 330.59 GHz were 22" and 23", respectively, which corresponds to 3300 and 3450 
AU at the distance to the target. The main beam efficiency was about 60 % around these 
frequencies. The system noise temperatures were typically 400 and 350 K for 12 CO and 
13 CO, respectively. For the backend, we used a 1024 channel digital auto-correlator, which 
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has a band width of 128 MHz and a resolution of 125 kHz, corresponding to 111 and 0.11 
km s -1 , respectively, at 345 GHz. The telescope pointing calibration was performed every 
1.5-2 hours by observing CRL 4211, and the resulting pointing accuracy was typically 3". 
For the absolute flux calibration, the integrated intensities of M1 7 SW and IRC 102 16 were 



measured two or three times per day and compared with those by IWang et al.l (11994I ). Data 
reduction and analysis were made with the NEWSTAR software package developed in NRO, 
which is a frontend of the Astronomical Image Processing System (AIPS) developed in the 
National Radio Astronomy Observatory (NRAO). 

We employed the five-point observation method in which spectra are obtained not only 
at a stellar position but also at four adjacent points 22" apart from the stellar position in 
the R.A. or Dec. direction. This method enables us to deduce the excess emission toward 
the central star by subtracting the average spectrum of the surrounding four spectra from 
the central spectrum, even if the star is located in extended components, such as an ambient 
cloud. 



3. Results of the Continuum Observations 

The 1.1 mm continuum map around V1094 Sco is shown in Figure [Th- Strong compact 
continuum emission was detected towards the stellar position. From the 2-d Gaussian fitting, 
the peak position of the emission is obtained to be 16 h 08 m 36 f 51, — 39°23'l'. / 6 in J2000, which 
agrees well with the 2MASS position of V1094 Sco within l'/5. The total flux density at 
1.1 mm was measured to be 272.1±7.6 mJy from the deconvolved cleaned map (see below), 
which is one of the strongest 1.1 mm emission among T Tauri stars in the Lupus molecular 
cloud. The negative contours around the emission are due to a point spread function (i.e., a 
beam; see Figure [lb) with a negative hole generated from the PCA cleaning method. 

Figure [lb shows a deconvolved map by subtracting the measured point spread function 
(PSF) from the emission via the clean algorithm. The PSF of the AzTEC map for V1094 Sco 
is shown in Figure [Tb. From the deconvolved map, we found that the emission is composed 
largely of a spatially unresolved component. The derived clean components are located 
around the star over two pixels (3" x a/2 ~ 4'.'2 in size) with its separation comparable to 
the pointing uncertainty estimated from the pointing measurements (Figure Hb). Therefore, 
we consider that the strong millimeter emission comes from a compact (r < 4'.'2/2 = 320 
AU) circumstellar dusty disk associated with V1094 Sco. This is s upported by the fact that, 



unlike the case of a stellar type flare seen in WTTS (lGudelll2002l ). V1094 Sco did not show 



any time variation of the flux density at millimeter wavelengths during the observing period. 



-6 - 



Figure [2] shows the SED of V1094 Sco. The flux density at 1.1 mm is unusually strong 
as the 1.1 mm t o infrared flux den s ity ratio is near unity and larger than those for T Tauri 
st ars in Taurus (IStrom et al.lll989t iBeckwith et al.lll990[ ). Adapting the extinction relation 
of iRieke fc Lebofskyl (119851 ). we measured the spectral index for the dereddened fluxes in 
the range from 2 to 10 /zm, a^io, to be —1.39 ± 0.10, where a = [dlog(AF\)/dlog(A)]. This 
value suggests that V1094 Sco should be categorized as a Class II object. The shape of the 
measured SED is unique in the sense that there are not distinct broadband spectral features, 
such as a large dip, from near-infrared (NIR) to millimeter wavelengths with a power-law like 
monotonic decrease. The index of the decrease is measured to be —0.98 ± 0.02 in the range 
from 2 to 1100 /xm. However, there is a small turnover around 10 /xm, i.e., the index at the 
longer wavelength, from 20 to 1100 jum, follows nearly the entire trend to be —0.93 ± 0.03 
while that from 2 to 10 fim is —1.39 ± 0.10. Furthermore, there is no dip feature in the NIR 
to MIR w avelength range cau sed by clearing material in the inner parts of the disk, such as 
GM Aur Jllughes et al.ll2009h . 



4. Results of the Molecular Line Observations 

The five-point spectra of the 12 CO and 13 CO J=3-2 emissions are shown in Figure [3j 
The significant 12 CO and 13 CO emissions were detected over the five-points around the star. 
There is a clear intensity gradient along the east-west direction in the 13 CO spectra as well 
as in 12 CO. In the 12 CO spectra, we could not find any symmetrical wing feature of an 
outflow with a velocity width of ~10 km s _1 which is often seen around a low-mass young 
stellar object. The 12 CO and 13 CO spectra are peaked around the LSR velocity of 4.9 km 
s -1 , whi ch is slightly shifted from the systemic velocity of the Lupus III main cloud (4.1 



km s 1 ; iTachihara et al.lll996l ). The peak intensity in Tmb scale and the velocity width in 
FWHM are 7.0 K and 0.66 km s -1 , respectively, on average over the five-points for 12 CO, 
and 0.7 K and 0.85 km s" 1 for 13 CO. 

From the intensity ratio of the 12 CO to 13 CO spectra, we estimate the optical depth 
of both lines toward the star. If we assume the same excitation temperature, the mean 
optical depths of 12 CO and 13 CO are obtained to be 16.7 and 0.19, respectively. Since the 
13 CO emission is optically thin in contrast to 12 CO, the associated gas with the star can be 
expected to be seen in 13 CO through the parent cloud as excess emission towards the star. 

By assuming the local thermodynamic equilibrium with the excitation temperature, T ex , 
the H2 mass of the gas in the beam solid angle toward the source, Mh 2 , can be estimated 
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from the intensity of the optically thin 13 CO J=3-2 line as follows: 
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M. 
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(1) 



where J(T ex ) = T /[exp(T /T ex ) — 1], T = hu/ks, Xrs co is the fractional abundance of 13 CO 
with respect to H 2 , D is the dis tance to the source , r is the optical depth, and tjmb is the 
main beam efficiency of ASTE ( iScoville et al.lll986l ). We adop t an excitation temper ature 
10 K which is the typical value in the ambient molecular cloud (jVilas-Boas et al.ll2000l ). and 
1.7 x 10 -6 for Xi3 CO by assuming that the abundance ratio of 12 CO/ 13 CO and the H 2 / 12 CO 
conversion factor are 60 and 10 4 , respectively. Our mass as well as column density are 
estimated to be 6.7 x 10 -3 M and 1.2 x 10 21 cm -2 by using equ ation (ED). The column 



density is much smaller than the averaged one in the Lupus 3 cloud ( iTachihara et al.lll996l ). 
suggesting that there is no large and massive envelope around V1094 Sco. 

To search for the compact excess emission towards the star, we made residual spectra 
of the 12 CO and 13 CO emission as shown in Figure HJ Each residual spectrum was obtained 
by subtracting the average of the surrounding four spectra from the central one. The excess 
emission above the 2a noise level over two velocity channels can be discerned in the 13 CO 
residual spectrum around the velocity of 6 km s _1 . The peak intensity and the velocity 
width in FWHM are measured to be 0.12 K and 1.43 km s -1 , respectively. Since the ve- 
l ocity of the excess em ission agrees well with the radial velocity of the star, 5.8±1.0 km s -1 
( IGuenther et al.ll2007l ). this may be the possible detection of the gas disk around V1094 Sco. 

The 12 CO residual emission was detected at Vlsr = 4.1 and 6.4 km s _1 whose peak 
intensities are 0.30 and 0.43 K in Tub, respectively. Although the 12 CO emission is shown to 
be optically thick, the double-peaked profile of the residual spectrum towards the star might 
suggest the presence of a rotating gas disk around VI 094 Sco. We tried to fit the rotating 
disk model ( iKitamura et al.lll993l ) to the both spectra to examine quantitatively. We could 
not, however, find any parameter set to reproduce the both residual spectra consistently. 
For example the 12 CO intensity can be reproduced for the larger rotating disk model. Such 
a disk model expects the stronger 13 CO intensity than the obtained and can not explain 
the both intensities compatibly. In addition, both velocity widths are not well fitted by the 
model simultaneously. Consequently, the residual spectrum in 12 CO is not likely to originate 
from the circumstellar gas disk. 
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Discussion 



5.1. Analysis of Spectral Energy Distribution 



The measured SED of V1094 Sco has the unique feature: there is no significant turnover 
and a monotonic decrease from NIR to millimete r. This feature is wel l descr ibed by a flat- 
tened disk for the dust component. According to IChiang &: Goldreichl ( 119971 ). the flattened 
disk has a lower temperature than that of a flared disk with the well-mixed dust because 
the central star can not illuminate the flattened disk efficiently. In this case, the SED is 
predicted to have an asymptotic power-law slope with an index of n = 4/3 from 30 //m to 
1 mm. The MIR to millimeter continuum radiation mainly comes from the optically thick 
interior of the disk. On the other hand, the SED of a flared disk rises in the MIR to far- 
infrared (FIR) regime owing to the strong emission from a hot surface layer and gradually 
falls down toward the millimeter regime, mainly tracing the cooler interior. 

To deduce disk parameters, we performed the least square fit ting to the SED of V10 94 Sco 
using a simple power-law disk model (same as model 1 case in lKitamura et al.l 120021 ). Our 
disk model is assumed to have a power-law form of the surface density and temperature 
variations over the entire area of the disk, S(r) and T(r), respectively: 



J 100AU 



and 



T(r) = T 1AV 



[100AU] 



[1AU] 



(2) 



(3) 



Here, r is the radial distance in units of AU, Siooau is the surface density at 100 AU, and 
T\A\j is the temperature at 1 AU. We applied the model to the data without the extinction 
correction because it was corrected in the model by c alculating the mas s opacity at each 
frequency from A v with the extinction curve (figure 1 of lAdams et al.lll988l ). The lower limit 
of the temperature is fi xed at 10 K which is the typical temperature of the Lupus 3 cloud 
( IVilas-Boas et al.ll2000f ). We adopt the effective temperature of the star of 4 255 K, the stellar 
mass of 0.78 M and the visual extinction toward the star of 1.75 mag ( jWichmann et al. 
19971 ). The outer radius of the disk, R on t, and the disk inclination angle, i, are fixed to be 
320 AU and 45 degree, respectively. The outer radius is derived from the pixel size of 3" 
as an upper limit because the emission was not spatially resolved. Stellar radius, i?*, inner 
radius of the disk, R in , Txau, Q and E 10 oau are the f ree parameters of the disk model. The 
model SEDs are computed and fitted to the o bservational data for the following four cases: 
p = 1 and 1.5 (minimum mass solar nebula by lHayashi et al.lll985l ) and /3 = and 1, where 



(3 is the power- law index of the dust mass opacity coefficient, k v = 0.1 x 



10 12 Hz^ 



cm 



-9 - 



( IBeckwith et al.lll990! ). In the fitting, we treated all the data points with the equal weights 
of the largest uncertainty because the SED fitting tends to be biased toward the data points 
with small uncertainties at shorter wavelengths. 

These model SEDs seem to agree with the observed SED as shown in Figure [2] and its 
reduced chi-squares are measured to be ~ 9. The derived parameters by the fitting are listed 
in Table [1] and [2j In Figure El we plotted the fitting results only for p = 1.5 case because the 
two best-fit curves with p = 1 and 1.5 are almost identical for each (3 value. Two remarkable 
disk features are derived in the fitting. First, the disk is very cold with a temperature of 99 
K at 1 AU. The fitted stellar luminosity of 0.85 L suggests that Tiau in the condition of 
radiative equilibrium is expected to be 269 K when the disk is in hydrostatic equilibrium and 
gas and dust are well mixed. Previous studies of the SED fitting ha ve found that most T 



Tauri stars in Taurus have a higher temperature of 100-400 K at 1 AU (IKitamura et al.ll2002 



Andrews fc Williams! 120051 ) . Our deri ved low temperature sugges ts a flattened dust disk as 
predicted by theoretical studies (e.g. IChiang fc Goldreichl 119971 ) . The second remarkable 
finding is that the disk is massive. Our estimated disk mass ranges from 0.03 to 0.12 M® 
and the ratio of the disk to stellar mass is 3-13%. These values are much higher than the 
typical values for T Tauri stars in the Taurus-Auriga region: median disk masses are 4 x 10~ 3 
and 2 x 10 -3 M m for CTTS and WTTS, respectively, and the median disk-to-star mass ratio 
is only 0.5 % (lAndrews fc Williamsll2005l ). It is yery rare to find a massive d isk in a quiescent 
accretion phase like a WTTS one. For example, Andrews fc Williams! (120051 ) have found that 
only 15 % of WTTSs are detectable at the submillimeter wavelengths. 

The above arguments are not affected much by either the disk inclination angle or 
the disk outer radius. If we change the disk inclination angle from 20 to 70 degrees with 
R out = 320 AU and p = 1.5, the temperature at 1 AU varies from 90 to 128 K, which is still 
low, and the index q remains around 0.68. The variation of the disk mass is less than 20 
%, remaining the massive disk. Furthermore, when we adopt a smaller R out than 320 AU, 
a higher surface density is required to reproduce the observed SED, resulting the massive 
disk. For example, in the i? out =110 AU case, Siooau becomes more than 5 times larger than 
that for the i? out =320 AU case and the disk mass increases up to 0.42 M , which remains 
massive. For a smaller opaque disk, we cannot accurately estimate the disk mass, though 
our interpretation of the massive disk is still valid. Therefore, we definitely conclude that 
the disk around V1094 Sco is colder than the usual flared disk and that the disk mass is 
much larger than the typical value of T Tauri stars in Taurus. 

In our SED fitting, we adopted a low threshold of 10 K in the temperature profile because 
the source is likely embedded in its parent cloud as shown by our 12 CO observations. This 
threshold causes an increment of flux densities around at 300 fim in the model SED because 
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the temperature of 10 K corresponds to 290 /im by Wien's law (see Figure |2]). To investigate 
the effect of the temperature threshold, we also analyzed the SED by using a disk model with 
a threshold of 2.7 K, the temperature for the cosmic background radiation. In this case, the 
observed SED can be well fitted only when the disk inclination angle is less than 30 degree. 
The fitted SED for p = 1.5 and (3 = 0.0 case is shown by a gray line in Figure [2] the best-fit 
curve can also explain the observed SED, and the disk mass becomes four times larger than 
that in the 10 K threshold case. Note that such a massive disk relative to the stellar mass 
is likely to be gravitationally unstable. On the other hand, Tiau and q decrease only by 2 % 
and 11 %, respectively. Consequently, the threshold in the temperature profile of the disk 
model does not change our interpretation of a cold massive disk around V1094 Sco. 



5.2. Rare Object to Provide the Initial Conditions of Planet Formation 



The cold and massive disk we found around V1094 Sco strongly suggests the presence of 
a flattened disk where a large number of small dust grains are settling. Since V1094 Sco seems 
to be in the quiescent phase of disk accretion, the majority of dust particles are expected 
to just settle into the midplane with dust growth through coalescence. Theoretically, the 
thin dust layer in the midplane fragments into numerous planetesimals by its gravitational 
instability, and as a result, a hole appears at the inner part of the disk to produce a dip 
feature at ~ 10 /im in SED. However, we could not find such a feature at the IR wavelengths 
in the SED of V1094 Sco, indicating that there is no inner hole created by the formation of 
planetesimals. Indeed, the strong continuum emission toward V1094 Sco that we detected 
suggests that a large number of small particles rather than larger planetesimals with sizes 
of ~ 10 km is the dominant component in the circumstellar disk. Consequently, the disk 
around V1094 Sco is most likely in the quiescent stage after the dust sedimentation in the 
midplane, but just prior to the planetesimal formation. 

V1094 Sco is a rare object in terms of the shape of its SED. To find a similar object with 
a massive cold disk to V 1094 Sco, we checked the available photo metric data at the MIR and 
millimeter wavelengths (IStrom et al.lll989t iBeckwith et al.lll990l ) and measured the flux ratio 
of MIR flux to millimeter flux. We found that only CY Tau, a Class II CTTS in Taurus, has 
a high intensity ratio of the 1 mm to 70 /im as well as to 25 /im (> 0.5), suggesting the strong 
millimeter flux despite its weak MIR flux and the SED with a monotonic slope from MIR to 
millimeter. The temperatu re at 1 AU and the di sk mass of CY Tau are estimated to be 87 K 
and 0.13 M , respectively (IKitamura et al.ll2002l ). We can roughly estimate the timescale for 
the duration of a cold and massive disk like those around V1094 Sco and CY Tau as follows. 
Given about 30 samples in Taurus whose flux densities at 1 mm and MIR were previously 
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determined, the timescale of the cold massive disk is estimated to be at least one order of 
magnitude shorter than the typical age of a T Tauri star (~ 1 x 10 6 ~ 7 yr), if we assume 
that new stars randomly form with a constant rate. Such a short timescale is consistent 
with the gravitational instability of the thin dust layer after the dust sedimentation at the 
midplane. It is suggested that the dust layer produced by settling to the equatorial plane 
bec omes gravitationally unstable and rapidly fragments into numerous planetesimal (<g;10 5 



yrs; 



Dominik et all 120071 ). 



It is very interesting that the physical properties of s uch a quiescent disk provide the 



initial conditions of the successive formation of planets (e.g. iKokubo fc Idall2002l ). In partic- 
ular, it is crucial to reveal the radial variation of the disk surface density with high-resolution 
and high-sensitive observations in future. For the dust coagulation, it is theoretically ex- 
pected that the index of the dust opac ity coefficient, 0, is lower than unity as a result of 
dust growth (IMiyake fc Nakagawalll993l ). We treated f3 as fixed parameter in the SED fitting 
because the photometric data at the millimeter and sub-millimeter wavelengths are limited. 
More (sub-)millimeter data are required to determine the value accurately. Additionally, 
high-resolution molecular line observations will reveal the disk kinematics and the initial 
conditions of the planet forming disk more quantitatively. 



6. Summary 

We have detected the strong 1.1 mm continuum emission toward the T Tauri star, 
VI 094 Sco, in the Lupus molecular cloud with AzTEC on ASTE. In addition, we have 
obtained the 12 CO and 13 CO J=3-2 spectra at the five positions around the source. The 
main results of our observations are summarized in the following. 

• Significant continuum emission was detected just toward the stellar position. The 
emission is not spatially resolved with the effective ASTE beam size of 35", indicating 
that a dust disk is associated with the star. The total flux density is 272 mJy at 1.1 
mm which is the strongest among T Tauri stars in the Lupus molecular cloud. 

• The 12 CO and 13 CO emissions were detected over the five-points around the star. The 
five-points spectra show no signature of a large-scale molecular outflow or a massive 
envelope. We made residual spectra by subtracting the average spectrum of the sur- 
rounding four spectra from the central spectra. The 12 CO excess emission was found 
and the 13 CO one was marginally discerned towards the star. The 13 CO excess may 
come from the gas disk while the 12 CO excess emission detected can not be explained 
by the gas disk. 
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• The observed SED shows two prominent features. First, there is no significant turnover 
with a constant slope from NIR to millimeter which is well reproduced by a flattened 
disk. Second, there is no signature of a dip at around 10 /im suggesting no inner hole 
in the disk. 

• We fit a power-law disk model to the SED when p is fixed to be 1.0 and 1.5, and (5 
is 0.0 and 1.0. The best fit parameters of the disk are Siooau =(1-25-4.97) g cm -2 , 
Tiau =99 K, q = 0.68, and M disk =(0.03-0.12) M . The disk temperature at 1 AU 
is lower than that of a flared disk in hydrostatic equilibrium. The disk mass is much 
higher than the typical values for T Tauri stars. 

• The cold and massive disk around V1094 Sco can be well interpreted as a partially 
settled disk of a large number of small dust grains. Together with no signature of an 
inner hole in the SED, it is most likely that the disk is in the quiescent phase after the 
dust sedimentation in the midplane, but just prior to the planetesimal formation. The 
V1094 Sco may be a rare sample to investigate the initial conditions for the formation 
of a planetary system. 
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Fig. 1. — (a) 1.1 mm continuum map of VI 094 Sco obtained with AzTEC on ASTE. The 
contour lines start at ±5<r with intervals of 5a where la is 5 mJy beam -1 . The negative levels 
are indicated by the dashed lines. The position of V1094 Sco is shown by the central cross, 
(b) The residual contour map by deconvolving the PSF via the clean algorithm is superposed 
on the continuum map (a) in grey scale. The contour levels and the grey scale are the same 
as in Figure Ufa). The shaded area indicates the positions of the clean components, (c) 
Point Spread Function (PSF) of the map (a) obtained with AzTEC on ASTE. The contour 
levels are ±10xn % of the peak value (n=l,2,3,. . .). The negative levels are indicated by the 
dashed lines. 



Table 1. Adopted parameters for SED fitting 



parameter 


value 


T eff [K] 


4255 


M* [M ] 


0.78 


A v [mag] 


1.75 


Rom [AU] 


320 


% [degree] 


45 


P 


1.0 & 1.5 





0.0 & 1.0 
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Frequency [Hz] 




10 100 
Wavelength [jum] 



1000 



Fig. 2. — Spectral Energy Distribution (SED) of V1094 Sco. The result of our observations is 
shown in the filled circl e. The filled squares in dicate the SED compile d from previous studies; 
the N OMAD catalog (IZacharias et al.l I2005T) . the DENIS d atabase (IThe Denis Consortium 
2005). the 2MASS point source catalog ( iCutri et al.l 120031 ) the IRAS faint so urce catalog 



fjMoshir et al.lll992l ). and the spitzer photometry (3.6, 4.5, 5.8, 8.0, and 24 fim ICieza et al. 



20071 ). The flux density at 70 /zm with Spitzer MIPS was measured by analyzing the archived 
data using the program MOPEX/APEX software packages. The IRAS flux densities at 12, 
60, and 100 /zm are upper limits and shown by the open triangles. The dashed lines show 
the stellar contributions and the solid line show the best-fitted curves in the case of p = 1.5. 
The red and blue lines show the calculations for /3 =0 and 1, respectively. A grey line 
indicates the fitting results for the model disk with a low threshold temperature of 2.7 K in 
the temperature profile when p = 1.5 and (3 = 0. 
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Fig. 3.— 12 CO (black) and 13 CO (red) J=3-2 five-point profile maps of V1094 Sco. The 
observed points are located at the stellar position, to the 22" north, south, east, and west 
from the center. The spectra are smoothed in resolutions of 0.3 and 0.5 km s _1 for 12 CO 
and 13 CO, respectively. The sha ded areas show the radial velocity of the central star with 
their error (iGuenther et al.l 120071 ). 
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Fig. 4. — Residual spectra of 12 CO and 13 CO towards V1094 Sco. The vert ical dotted lines 
and t he shaded areas show the radial velocity of the star and its uncertainty ( iGuenther et al. 
2007h . 



Table 2. Results of SED fitting 



parameter 




value 


P 


0.0 1.0 


V 


1.0 


1.5 1.0 1.5 


R* [AU] 




0.008 






0.85 


Rin [AU] 




0.02 


Siooau [g cm" 2 ] 


1.25 


1.47 4.31 4.97 


T 1AU [K] 


99.0 


99.0 99.3 99.4 


Q 




0.68 


M disk [M Q ] 


0.03 


0.04 0.10 0.12 
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